Diffusion regimes in Levy nights with trapping 
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The diffusion of a walk in the presence of traps is investigated. Different diffusion regimes are 
obtained considering the magnitude of the fluctuations in waiting times and jump distances. A 
constant velocity during the jump motion is assumed to avoid the divergence of the mean squared 
displacement. Using the limit theorems of the theory of Levy stable distributions we have provided 
a characterization of the different diffusion regimes. 
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I. INTRODUCTION 

There are many physical situations where diffusion 
takes place under the presence of traps, trapping dif- 
fusion. Examples are found in electronic conduction in 
amorphous semiconductors and quasicrystals [fij, atomic 



, tracer diffusion in 
,|).The existence of 



diffusion in glass like materials 
living polymers J3|, and more 

traps is generally modeled through a probability density 
of waiting times between successive steps in the walk, 
continuous-time random walks (CTRW) ||. The theory 
of CTRW has been extensively developed in the litera- 
ture, using the generating function methods |g{ or simple 
statistical reasoning Q . The existence of a wide distribu- 
tion of waiting times leads to a subdiffusivc regime where 
the mean squared displacement grows slower than time. 
On the other hand, there are physical situations where 
the time between successive jumps may be considered 
constant, but the distribution of jump distances is wide 
gggg, (Levy flights) |. ' 

Levy flights are also observed in a large variety of phe- 
nomena, for instance in chaotic dynamical systems ftl\p \, 
turbulent PJTc|] flow and self-gravitating systems Jll|. 
The existence of a wide distribution of jump distances 
leads to a superdiffusive regime, where the mean squared 
displacement grows faster than time. Moreover, Levy 
flights generates fractal structures in space 0. Nev- 
ertheless, more complex behaviors are expected in sys- 
tems where both the distribution of waiting times and 
jump distances are wide. For instance, Schulz p2| stud- 
ied an anomalous Drude model for transport properties 
of quasicrystals. He assumed that the walks move with 
an anomalous speed v a ~ t a between collisions, and a 
power law p(r) ~ T - 1 -^ distribution of time between 
collisions. In this way, he obtained the phase diagram 
(/i, a) of the system, which is divided in different regions. 
While the anomalous velocity may mimic the existence of 
traps, it is introduced artificially and cannot be derived 
from a Hamiltonian. In [ p^| , Klafter et al introduced a 
stochastic description of anomalous transport phenom- 
ena and found different behavior of the mean square dis- 



placement of CTRW. Zumofen and Klafter g studied a 
one-dimensional map which exhibits intermittent chaotic 
behavior with coexisting laminar and localized phases, 
and analyzed them in terms of Levy statistics. 

A more consistent approach was presented by Klafter 
and Zumofen Ej . They studied the diffusion in a Hamil- 
tonian system in terms of the CTRW formulation, and 
considered wide distributions of waiting times and jump 
distances. However, they did not provide a complete 
characterization of the different diffusion regimes that 
may be obtained. In a recent letter Fogedby p5[ consid- 
ered Levy flights in the presence of a quenched isotropic 
random force field studying the interplay between the 
built in superdiffusive behavior of the Levy flights and 
the pinning effect of the random environment. 

In the present work we study the behavior of a random 
walk in the presence of traps. Disorder is introduced in an 
annealed way, through power law distributions of waiting 
times p(r) ~ T -i-°<™ and jump distances p{x) ~ x ~ a:c . 
We use simple statistical reasoning close to the approach 
developed by Bouchaud and Georges [H instead of that 
of the work by Klafter and Zumofen |3 , which use the 
CTRW formulation. The present formalism contains as a 
fundamental tool the theory of Levy stable distributions, 
introduced by Levy |]l6f and developed by other authors 

OS 



II. THE MODEL 

We consider a random walk on a lattice, such that 
the particle has to wait for a time t w on each site be- 
fore performing the next jump. The waiting time is a 
random variable independently chosen at each new jump 
according to a distribution p(r w ). We also assume that 
the waiting time is not correlated to the length of the 
jump x, which is distributed according to p(x) . The diffu- 
sion process will be characterized by the scattering func- 
tion F(k,i), the Fourier transform of the diffusion front. 
Other properties like the diffusion front and the mean 



1 



squared displacement can be derived from this function. 
For instance, the mean squared displacement is given by 



(x 2 (t)) = 



d 2 F 



dk 2 



(1) 



fe=0 



Let iV be the number of steps performed by a walker 
during time t. N is, in general, a random variable which 
depends on the duration of the jumps and waiting times. 
The scattering function can thus be expressed as a sum 
over all possible jumps during time t 



F(k,t)= dNF(k,N)P(N,t), 



(2) 



where F(k, t) is the scattering function of the same prob- 
lem, but considering regular duration of the jumps and 
no waiting time and P(N, t) stands for the probability 
distribution of N jumps at a fixed time t. 



A. Mean squared displacement after TV steps 



The total displacement after N steps is given by 



(3) 



3 



2 







- Vc 


1 1 lb 


I ■ 


Vb 


Ilia / 


II 








IV 





1 



2 



3 



a 



w 

FIG. 1. The phase diagram (a w , a x ). The behavior in the 
different regions is as follows: I.- normal diffusion; II.- LTT 
with exponent a x and superdiffusion; III.- LTT with expo- 
nent a m a) superdiffusive, b) normal diffusion; IV.- LTT with 
exponent a x and ballistic motion; V.- LTT with exponent 
a w , a) superdiffusive, b) and c) subdiffusive. See text for a 
detailed description. 



In the right hand side we have a sum of mutually inde- 
pendent random variables with the common distribution 
p(x), with zero mean. The limit distribution for large N 
will be a stable Levy distribution 0|L8fl, i.e. 



X N =l*N 1/a * 



(4) 



where = denotes that random variables in both sides have 
the same distribution, I* is a characteristic value and u 
follows the symmetric Levy distribution L ax Q(u). The 
canonical (Fourier transform) representation of Levy sta- 
ble laws is (for a ^ 1) 



FT[L aP ](k) =exp 



\k\ a 1 



k it 
t— rip tana— 
\k\ 2 



(5) 



where a and (3 are real numbers defined in the intervals 
< a < 2 and -1 < < 1. The case a = 2 and (3 = 
corresponds with the Gaussian distribution, which de- 
cays faster than any power law for large arguments. On 
the contrary, all Levy distributions, except the Gaussian, 
have the asymptotic behavior for u^> 1 p7 18 



Then, from eqs. (Q) and (||) it follows that 

F(k,N) = exp[-(|fc|P) aa =7V], 



exp[ 

If p(x) has finite variance a then I* — a and a 
If p(x) « l^lxl- 1 -^, with < fi < 2 then I* ~ I 
a x = [i. 



(6) 



(7) 

= 2. 

o and 



B. Number of steps after time t 

On the other hand, the number of steps after time t is 
given by 



N 



N 



(8) 



where T X i are the duration of the jumps. If we assume 
that during the jump motion the walker moves continu- 
ously at a constant velocity v and changes directions at 
random then t x = i> -1 |x| and p(r x ) = 2vp(\x\/v). In the 
right hand side of eq. (0) we have two sums of indepen- 
dent random variables, with common distribution p(r w ) 
and p(r x ), respectively. The limit distributions for large 
N will follow Levy distributions |l7|Jl§f|, i.e. 



t = tN + T*N 1/a ™ Ul + T* x N 1/a > 



1t 2 . 



(9) 



where u\ and u-i follows the Levy distribution i a „i(iti) 
and Lq, x 1(1x2), respectively. The first term in the right 
hand side appears only if p(r w ) or p(r x ) have finite mean, 
and t is given by the sum of the finite means. If p(r w ) 
has finite variance a then r* = a and a w — 2, while 



t^t- 1 -^ (0 < a* < 2) then r* 



To and 
= a and 



if p(r w ) 

a w = [i. If p(t x ) has finite variance a then r* 
a x = 2, while if p(r x ) « Tq x tz}-' 1 (0 < /i < 2) then 
t* ~ t and a x = ^.From eqs. (g), (0) and (Jsj) it follows 
that 
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F(k,t) = J I du 1 du 2 L awl (u 1 )L axl (u 2 )ex.p[-{kl*) a 'N]. 

(10) 

where the functional dependence of N with m, U2 and t 
is determined from eq. (||) . Next we analyze the behavior 
of the scattering function, defined trough this expression, 
for different values of a x and a w . With this purpose we 
have divided the corresponding phase diagram in five re- 
gions, as it is illustrated in fig. 1. 



III. RESULTS 



In region I both the distribution of waiting times and 
jump distances have finite variance. Hence, the random 
variables ui and 112 will follow a Gaussian distribution. 
For large N, we can therefore neglect the last two terms 
in the right hand side of eq. (||), obtaining N « tjr. 
Then eq. ([l0|) reduces to 



F(k,t) = exp[-{\k\l*)H/r]. 
Moreover, using eq. (Q) one obtains 

~t. 



(11) 



(12) 



We thus find normal or classical diffusion in this region: 
the scattering function decays exponentially with time 
and the mean squared displacement grows proportional 
to time. 

In region IV the distribution of jump distances is quite 
wide, and even wider than the distribution of waiting 
times. Hence, it is expected that the third term in 
eq. (^) gives the major contribution, thus obtaining 
N » (i/r*) Q -u^ Qx . Substituting this result in eq. $sfy, 
and expanding the exponential inside the integral, one 
obtains 



F(k,t) 



^-^ r(l + not 

n=0 y 



("!)" ( l)na. 

Tk 



(13) 



where Tf. — T*/\k\l*. This series has an infinite radius 
of convergence and, therefore, can be taken as a se- 
ries expansion for the scattering function in this region. 
For small times (t <C Tk) the scattering function follows 
a stretched exponential decay, with stretched exponent 
a x . On the contrary, for ( > the relaxation becomes 
slower than an exponential. Taking L a ^ i{u\) ~ 
from cq. @ it follows that F(k,t) ~ t- a * . Hence, for 
long times the scattering function follows a long time tail 
(LTT), with an exponent a x smaller than one. In this 
case, the mean squared displacement determined from 
eqs. (Q) and ( |l3| ) is not finite. This is a consequence of 
the divergence of the second moment of the distribution 
of jump distances. Nevertheless, the total displacement 
at time t cannot be larger than vt and, therefore, there is 



a cutoff k c ~ t _1 for small values of k. Thus, to avoid the 
divergence of the mean squared displacement we evaluate 
cq. (|jj) in k — k c instead of k = 0. In this way we obtain 



(x 2 (t)) ~t 2 , in IV. 



(14) 



The motion of the walk is in this case of ballistic type. 
The behavior in this region has been investigated by dif- 
ferent authors, which consider random walk motion due 



to a periodic potential |14|. 

In region V we can use a similar reasoning as in region 
IV, but in this case the dominant term will be the second 
one, associated to huge fluctuations in the waiting time. 
In this way we obtain 



00 



n— 



(15) 



where = T^{\k\l*)~ a:c / am . This expression is similar 
to eq. ([l3|). Besides, the asymptotic behaviors for short 
and long times are the same, replacing a x by a w . The 
main difference is observed in the /c-dependence of the 
relaxation time Tk, which is manifested in the temporal 
dependence of the mean squared displacement, 



( X *(t)} 



2— a x +o:„ 



in Va and Vb; 
in Vc. 



(16) 



The difference between region Va and Vb is that, in the 
former the mean squared displacement grows faster than 
time and, therefore, the system is in a superdiffusivc 
regime, while in the second one it grows slower than time 
and the system is in a subdiffusive regime. In region Vc 
there is also a subdiffusive regime, which has been investi- 
gated by different authors, using the generating function 
formulation Q or simple statistical reasoning How- 
ever, in region Vc a x > 2 and, therefore, the spatial tra- 
jectories of the walk will not be self-similar as in regions 
Va and Vb. 

In region II and III the fluctuations in jump distances 
and waiting times, respectively, are still Levy type, but 
are no so strong as in regions IV and V, respectively. In 
these cases one cannot neglect the first term in the right 
hand side of eq. (||) and, hence, there is some reminiscent 
of normal diffusion behavior. 

In region II the distribution of jump distances is wider 
than the distribution of waiting times, and both have fi- 
nite mean. We thus expect that the third term in the 
right hand side of cq. (O) gives the major contribution 
to the fluctuations in N, and the second one may be ne- 
glected. Even with this simplification we cannot solve 
eq. (^|) analytically, however the following asymptotic 
behaviors are obtained 



N : 



t/r u <C u c (t), 

{t/T*) a * u->u c {t); 



(17) 



where u c (t) ~ (t/r)" 1 . For small times, substituting eq. 
( |j"7| ) in cq. ( JTo| ) one obtains 



3 



F(k,t) = exp[-(kl*) a *t/T}. 



(18) 



The relaxation for small times is the exponential like in 
the normal diffusion case, region I. However the mean 
squared displacement does not grow linearly with time, 



(* 2 (*)>~t 3 



in II. 



(19) 



but is characteristic of a superdiffusion regime. For long 
times, from eqs. ( |l0| ) and eq. (|l7|), and the asymptotic 
expansion for large arguments of Levy distributions in 
eq. (H), it follows that F(k,t) ~ t~ ai . For long times 
the scattering function decays slower than an exponen- 
tial following a power tail, like in region IV, but with a 
larger exponent. 

In region III we can use a similar reasoning as in region 
II, but now neglecting the fluctuations in the jump dis- 
tances in relation to the fluctuations in the waiting times. 
Again, even with this simplification we cannot solve eq. 
(^) analytically, however the following asymptotic behav- 
iors are obtained 



N 



t/r 
(t/r* 



u < u c (t), 
u > u c (t); 



(20) 



where now u c (t) ~ (t/r) Q ™ . For small times, substituting 
eq. (|2^) in eq. (|i"o| ) one obtains an exponential relaxation 
as in eq. (|l8|). However, the characteristic exponent a x 
can be now larger than two leading to different behaviors 
for the mean squared displacement, 



(x 2 (t)) 



in Ilia; 
in 111b. 



(21) 



Hence, in subregions Ilia and Illb we found superdiffu- 
sion and normal diffusion behavior . respectively. For long 
times, from eqs. ( |Io|) and eq. (|2C|), and the asymptotic 
expansion for large arguments of Levy distributions in 
eq. (^|), it follows that F(k,t) ~ t~ Uu > . For long times 
the scattering function decays slower than an exponential 
following a power tail, like in region V, but with a larger 
exponent. 



IV. SUMMARY AN CONCLUSIONS 

In summary we have investigated the diffusion behav- 
ior in the whole plane (a w , ot x ), which has been divided 
in five regions considering the magnitude of the fluctu- 
ations in waiting times and jump distances. A constant 
velocity during the jump motion was assumed to avoid 
the divergence of the mean squared displacement. Using 
as a fundamental tool the limit theorems of the theory of 
Levy distributions we have provided a characterization 
of the different diffusion regimes. 

We conclude that the diffusion behavior cannot be just 
classified as normal diffusion, superdiffusion and subdif- 
fusion. This classification only takes into account the 
temporal dependence of the scattering function for short 
times, while the long time decay may be different. For 



instance, in regions II and Ilia we have similar superdif- 
fusive regimes, but the long time behavior is determined 
by the distribution of jump distances, in the first case, 
and by the distribution of waiting times in the second 
one. 
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